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Abstract 

In this paper we present a study of the effects of two different single-ion anisotropies in the phase diagram and in the 
compensation temperature of the mixed spin-2 and spin-5/2 Ising ferrimagnetic system. We employed the mean-field theory 
based on the Bogoliubov inequality for the Gibbs free energy. Also we found the Landau expansion of the free energy in 
the order parameter to describe the phase diagrams. In the plane critical temperature versus single-ion anisotropic the phase 
diagram displays tricritical behavior. The critical and compesation temperatures increase when the single-ion anisotropies 
increase. 
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I. INTRODUCTION 



In the last five decades the Ising model has been one of the most largely used to describe critical behavior of 
several systems in the nature. In particular, in the physics of the condensed matter it is important to describe critical 
behavior and other thermodynamics properties of a variety of physic systems (disordered system, spins glass, random 
field Ising model, etc.). Recently, several extensions have been made in the spin— 1/2 Ising model to describe a wide 
variety of physic systems. For example, the models consisting of mixed spins of different magnitudes are interesting 
extensions, which are the so-called mixed-spin Ising models. 

The researches in ferrimagnetic materials are of great interest due to their possible technological applications, as well 
as of the academic point of view. These materials are modeled by mixed-spin Ising model, which can be constituted 
by several combinations of spins (spin-1/2, spin-1), (spin-1/2, spin-3/2), (spin-1, spin-3/2), (spin-2, spin-3/2), (spin-2, 
spin-5/2). The interest in studying magnetic properties of some types of ferrimagnetism, namely the molecular- based 
magnetic materials f^^W\, is due to its less translational symmetry than to their single-spin counterparts since they 
consist of two interpenetrating sublattices. The bimetallic chain complex MnNi(EDTA) — 6H2O is an example of an 
experimental mixed-spin system 

Additionally, there are many studies of mixed-spin Ising systems for explanation of the physical properties of 
disordered systems. This theme has been a great challenge in statistical mechanics. In this sense, in the last years 
there has been great interest in the study of magnetic properties of systems formed by two sublattices with different 
spins and with different crystal field interactions|t5i]. 

One of the earliest and simplest of these models to be studied was the mixed-spin Ising system consisting of spin-1/2 
and spin-S' {S > 1/2) in an uniaxial crystal field 0, 0|- From the theoretical point of view, such systems have been 
widely studied by a variety of approaches, for instance, effective-field theory ^Snlfi], mean-field approximation [l7j . 
renormalization-group technique [l^ and Monte-Carlo simulation (20l - [28l | and exact solutions for the mixed spin-1 
and spin-S Ising model in an uniaxial crystal field (29l. |30|. 

Recently, the goal is to extend the investigations to a more general mixed-spin Ising model with one constituent 
spin-1 and, in the simplest case, the other constituent spin-3/2. Abubrig et al. [31] presented a study of mean-field 
theory to determine the effects in the phase diagram of different crystal fields of the mixed spin-1 and spin-3/2 Ising 
system. They also showed some outstanding features in the temperature dependence of the total and sublattice 
magnetization. 

In this paper, we studied the effects of two different single-ion anisotropics in the phase diagram and in the 
compensation temperature of the mixed spin-2 and spin-5/2 Ising ferrimagnetic system. The outline of the remainder 
of the paper is as follows: In Section II, the model is introduced and we have obtained analytical expressions for 
free-energy, equations of state. We also have found the Landau expansion of the free energy in the order parameter. 
In Section III, we present the results and discussions about the phase diagrams and the compensation temperature. 
Finally, in Section IV we present our conclusions. 

II. THE MODEL AND CALCULATION 

The mixed-spin ferrimagnetic Ising system consists of two interpenetrating square sublattices (A and B) with spin 
S'"^ = 0, ±1,±2 and spin = ±1/2, ±3/2, ±5/2. In each site of the lattice there is a single-ion anisotropic {Da 
in the sublattice A and Db in the sublattice B) acting in the spin-2 and spin-5/2. The system is described by the 
following Hamiltonian model: 



where the first term represents interaction between the nearest neighbors spins in the sites i and j located in the 
different sublattices A and B. J \s the magnitude of this interaction, and the sum is over all nearest neighboring pairs 
of spins. The second and third terms represent the single-ion anisotropics at all points of the sublattices A and S, 
respectively. The sums are performed on A^/2 spins of each sublattice. 

In order to derive the analytical expressions for free-energy and equations of state, we employed the variational 
method based on the Bogoliubov inequality for the Gibbs free energy 



where G{H) is the free energy of H, and Go('Ho) is the free energy of a trial Hamiltonian 'Ho{r]) depending on 
variational parameters. (■ • ■ )o denotes a thermal average over the ensemble defined by T-Lo{r]). To facilitate the 
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calculations, we choose the simplest trial Hamiltonian, which is given by 



(3) 



where rjA and ijb are variational parameters related to the two different spins. Through this approach, we found the 
free energy and the equations of state (sublattice magnetization per site ttia and ttib)' 
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where g — , P = ^/kBT, N is the total number of sites of the lattice and z is the coordination number. The 

M\„ ^^A ^„ — /qB\ 



sublattice magnetization per site toa and ms are defined by niA = ('5'j )o and tub = {Sj)o, thus 

2 sinh (2/3r/A) + exp (-3/31?^) sinh {(Stja) 
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Minimizing the free energy in terms of the variational parameters rjA and rjB-, we have obtained 

r/A = JziriB , ?7_B = JziTiA- 



mB = \ 
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Therefore, we can get the generic mean-field equations (HUT]), which provide the magnetic properties of ferrimanetic 
system. Since these equations ([5][7]) have in general several solutions for the pair (m^ and ms), the stable phase will 
be one which minimizes the free energy. The detailed phase diagram is determined by numerical analysis, but some 
features of the phases diagram can be obtained analytically. Therefore, close to the second-order phase transition 
from an ordered state (toa ^ and ms ^ 0) to a disordered state (toa = and = 0). The magnetization mA and 
niB are very small, so we can expand the equations (j4]|6)) to obtain a Landau-like expansion: 



where the expansion coefficients are given by 
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III. RESULTS AND DISCUSSIONS 



The phase diagrams were constructed according to the fohowing routine: i) numerical solutions of A2 = and 
A4 > provides second-order transition lines, ii) A2 = 0, A4 = and Aq > determines the tricritical points, 
iii) The first-order transition lines are determined by comparing the corresponding Gibbs free energies of the various 
solutions of equations ([5]) and ^ for the pair (m^, m^). Even so, we have also analysed that > in all T, Da, Db 
space. 

The particular case Da = Db ~ 0, the critical temperature is determined by taking vtla — > and nis ^0, or 
= and Da^Db = 0, thus fcsTc/ J = 9.6609. 



A. The ground-state 

The ground-state phase diagram (see Fig. 1) is determined from the Hamiltonian ^ by comparing the ground-state 
energies of the different phases. At zero temperature, the structure of the ground state of the system consists of six 
phases with different values of {mA,TnB, Qa, qs}, namely the ordered ferrimagnetic phases 
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Figure 1: Ground-state phase diagram of tiie mixed spin-2 and spin-5/2 Ising ferrimagnetic system with two different single- 
ion anisotropies Da/\J\ and Db/\J\- The six ordered phases are represents by Oi{— 2, 5/2, 4, 25/4}, 02{ — l, 5/2, 1, 25/4}, 
03{-2, 3/2, 4, 9/4}, 04{-l, 3/2, 4, 9/4}, 05{-2, 1/2, 4, 1/4}, Oej-l, 1/2, 1, 1/4}. 



B. Phase diagram 

The phase diagrams are analysed in the {Da/ \ J\ , kBT/\J\) and {Db/ \ J\ , fc^T/l J|) planes, where the numerical 
solutions of ^2 = with A4 > provides second-order transition lines. The tricritical points are determined by 
numerical solutions of A2 = 0, A4 = Q with Aq > 0. The first-order transition lines were determined by comparing 
the corresponding Gibbs free energies of the various solutions of equations ([5]) and ([6]) for the pair (ttia, to_b). 

In Fig. 2, we display the phase diagram of kBTc/\J\ versus Da/\J\ for selected values of Db/\J\- In regions of 
high temperatures, for all positives and negatives values Da/\J\, and for any values of Db/\J\, the phase diagram 
shows only second-order phase transitions, which are indicated by solid lines. For values of Db/\J\ > 10.0, all 
second-order lines end in the same tricritical point given by {kBTt/\J\ — 4.5254, _D^/|J| = —9.3200). However, for 
values of /\J\ < —10.0, now in low temperatures, also all second-order lines end in same tricritical point given by 
{kBTc/\J\ = 0.8972, Da/\J\ = —1.8667). One heavy dotted curve connects these two tricritical points {kBTc/\J\ — 
4.5254, Da/I J| = -9.3200) and (fcsTc/lJl = 0.89716,1)^/1^1 = -1.8667). This curve separates the region with 
second-order phase transition of the region with first-order phase transition. In the region of low temperatures and 
for all values of Da/\J\ and Db/\J\, the phase transitions are of first-order. Thus, in this space {Da/ \ J\ , fesT/j J|), 
the system presents tricritical behavior. 

Additionally, the diagram shows that when Db/ \ J\ +00, the mixed spin Ising system behaves like a two-level 
system since the spin-5/2 behaves like = ±5/2. Nevertheless, in case that Db/ \ J\ — —00, the = ±5/2 states 
are suppressed and the system becomes equivalent to a mixed spin-1/2 and spin-2 Ising model. So, this is the reason 
that the coordinates of the tricritical point in the limit of large positive Db/ \J\ are well higher than those for large 
negative Db/ \J\- For the special case with equal anisotropic fields {Da/ \ J\ ^ Db / \ J\ = 0), the critical temperature 
is kBTc/\J\ = 9.6609, and the location of the tricritical point is kBTt/\J\ = 3.3833, Da/I^I = -6.2150, for z = 4. 

In Fig. 3, it is shown the phase diagram of kBTc/\J\ versus Db/\J\ for various values of Da/\J\- In the case of 
Da/\J\ > —1.8842 the phase transitions are only of second-order (solid lines) for any values of Db/\J\- The value of 
the critical temperature increases when Db/\J\ and Da/\J\ also increase. Still, in range —1.8842 < Da/\J\ < —9.0, 
the phase transitions are of second-order in the high temperatures region (solid lines) and of first-order (light dashed 
lines) in the low temperatures region. One heavy dotted curve of tricritical points separates the second from the 
first-order transition lines. 

One interesting feature shown in the diagram of the Fig. 3, refers to the fact that the phase transitions are only of 
first-order when the values of Da/\J\ < —9.0. All lines that start below the heavy dotted line will necessarily be only 
of first-order phase transition. Again, in this space (Dr/I J|, A:bT/| J|), the system presents tricritical behavior. 

These results may be compared to those obtained in the paper [l3|, which it has an error in equation (4) and 
consequently in the equation (8). These errors led to very different phase diagrams at temperatures above zero. 
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Figure 2: Phase diagram in the (Da/| Jj, fcsT/| J|) plane for the mixed-spin Ising ferrimagnet with the coordination number 
2 = 4, for several values of Db/\J\- The solid and light dashed lines, respectively, indicate second- and first-order phase 
transitions, while the heavy dashed line represents the positions of tricritical points. 




Figure 3: Phase diagram in the (Ds/j J|, fc_gr/| J|) plane for the mixed-spin Ising ferrimagnet with the coordenation number 
2 = 4, for several values of Da/\J\- The solid and light dashed lines, respectively, indicate second- and first-order phase 
transitions, while the heavy dashed line represents the positions of tricritical points. 



C. Compensation temperature 

The present mixed-spin system can exhibit compensation points, and to show this fact, we will consider J < 0. 
The signs of the sublattice magnetizations are different since we are taking into account that in the ferrimagnetic 
case the system consists of two interpenetrating square sublattices (A and B) with spin-2 and spin-5/2. Thus, it is 
possible that there are a compensation temperature Tcomp with Tcomp < 7c. So, the total magnetization per site 
M = {mA + rnB)l2 is equal to zero, same that mA ^ and txib ^ 0. 

In Fig. 4, we present the diagram Tc and Tcomp versus Da/\J\, and for some selected values of Db/\J\- The 
diagram shows that there are compensation points in the range 2.5339 < Db/\J\ < —10.0 and for Da/ \J\ > —2.50, 
which are indicated by dotted lines, while solid lines indicate critical temperature. The inset in Fig. 4 exhibits some 
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Figure 4: The critical Tc and compensation Tcomp temperatures as a function of the single-ion anisotropy Da/\J\, and for 
several values of Db/\J\- The solid and dotted curves represent the critical temperature and compensation temperature, 
respectively. The inset shows a magnification of the region closed where the compensation and critical temperatures are 
together for Db/ \J\ = —10.0. Temperatures are measured in units of \ J\ /ks- 

compensation points in case of Db/ \ J\ = —10.0. The Fig. 5 exhibits the diagram Tc and Tcomp versus Db/\J\, for the 
some selected values of Da/\J\- This figure only confirms the information shown in Fig. 4, ie, there is compensation 
temperature in the range 2.5339 < Db/\J\ < —10.0 and Da/ \ J\ > —2.5. 
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Figure 5: The critical Tc and compensation Tcomp temperatures as a function of the single-ion anisotropy Db /\J\, and for several 
values of Da/\J\- The solid and dotted curves represent the critical temperature and compensation temperature, respectively. 
Temperatures are measured in units of |J| /fc^. 



IV. CONCLUSIONS 

In this paper, we have studied the effects of two different anisotropics in the phase diagram and in the compensation 
temperature of the mixed spin-2 and spin-5/2 ferrimagnetic Ising system by using the mean field theory based on 
the Bogoliubov inequality. The phase diagrams are shown in the critical temperature versus single ions anisotropics 
plane. The system presents tricritical behavior, ie, the second-order phase transition line is separated of the first-order 
transition line by a tricritical point. Additionally, we also observed compensation temperatures. So, in conclusion, we 
can say that such a system may exhibit tricritical behavior and compensation temperatures due to the two different 
anisotropics. 
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